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We consider multiclass feedforward queueing networks with first 
in first out and priority service disciplines at the nodes, and class 
dependent deterministic routing between nodes. The random behav- 
ior of the network is constructed from cumulative arrival and service 
time processes which are assumed to satisfy an appropriate sample 
path large deviation principle. We establish logarithmic asymptotics 
of large deviations for waiting time, idle time, queue length, depar- 
ture and sojourn-time processes in critical loading. This transfers 
similar results from Puhalskii about single class queueing networks 
with feedback to multiclass feedforward queueing networks, and com- 
plements diffusion approximation results from Peterson. An example 
with renewal inter arrival and service time processes yields the rate 
function of a reflected Brownian motion. The model directly captures 
stationary situations. 

1. Introduction. A frequently considered asymptotic for analyzing the 
complex behavior of queueing networks deals with the situation when the 
queues approach a critical load: Based on a functional central limit the- 
orem for the driving processes, the weak convergence of suitably normal- 
ized network behavior to reflected (fractional) Brownian motion can be 
proved [13, 17, 21, 29, 34]. Here we speak of a reflected process when it 
is restricted to the positive orthant through a Skorokhod map [9, 14]; see 
Remark 5 after Theorem 2 for a definition. This heavy traffic convergence 
has been one motivation for a large number of publications about large 
deviations of reflected Gaussian processes [1, 20, 23, 25, 27, 28]. However, 
the weak convergence of queueing processes to a reflected Gaussian process. 
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is not appropriate for justifying a relation between large deviations of the 
queueing processes and large deviations of the reflected Gaussian limit pro- 
cess: Weak convergence focuses on expectations, whereas large deviations 
characterize logarithmic tail asymptotics. In this work we partially fill this 
gap by providing a framework for large deviation asymptotics of multiclass 
feedforward queueing networks with first come first serve (fifo) and prior- 
ity service disciplines at the nodes in critical loading. Feedforward means 
that the nodes can be numbered in such a way that the sequences of sta- 
tions which are visited by customers have strictly increasing numbers. We 
consider deterministic class dependent routing of customers between nodes. 

In the setting of single class queueing networks such a framework has 
been developed by Puhalskii [32]. There, a convergence toward critical load- 
ing is combined with a moderate deviation principle for the centered input 
processes. In this work a moderate deviation principle is regarded as a large 
deviation principle with a scaling close to a functional central limit theorem 
such that the rate function for the tail probabilities of a Gaussian process 
appears [6, 26, 35]. In this way Puhalskii obtains the characterization of 
logarithmic tail asymptotics for the appropriately scaled behavior of single 
class queueing networks when the queues approach a critical load. The rate 
function in the examples of Puhalskii [32] also appears in large deviation 
principles for reflected Brownian motion. Puhalskii's approach and its car- 
ryover to multiclass feedforward networks, presented here, is able to take 
advantage of both the simplified heavy traffic model captured through a 
Skorokhod map and the explicit rate function of a Gaussian process. 

Puhalskii's work addresses random routing and feedback between the 
queues which is not permitted in our work where the feedforward struc- 
ture allows for an iterative propagation of certain properties through the 
queues of the network. In our situation the difficulties come rather from 
the multiclass fifo service discipline. We note that in technical applications 
(e.g., the analysis of broadband communication networks) one is frequently 
faced with multiclass queueing networks where (at least in large parts of 
the network) feedback is not relevant, but the mixture of fifo and priority 
service disciplines. 

The starting point for our results is a large deviation principle for the se- 
quence of centered input processes which characterize the exogenous arrival 
and service times at the nodes. We then make use of the model for the be- 
havior of multiclass feedforward queueing networks developed in [21]. There, 
it has been shown that close to a homogeneous critical load situation, this 
behavior can be approximated through a multi-dimensional Skorokhod map. 
We use this pathwise heavy traffic convergence to deduce sample path large 
deviation principles for the sequences of idle-time, workload, queue length, 
departure and sojourn time processes of the queueing network induced by 



LARGE DEVIATIONS FOR FEEDFORWARD QUEUEING NETWORKS 



3 



the sequence of input processes. Besides the mentioned pathwise heavy traf- 
fic convergence, the proof of this result reUes on a version of the contraction 
principle provided in [31]. 

In order to obtain these results, the mean arrival and service rates must 
be chosen in such a way that a critical load is approached at a certain speed 
at every node. In addition, our heavy traffic condition requires that the 
mean rates in the sequence of input processes grow to infinity. A moderate 
deviation principle can provide the basis for satisfying this heavy traffic con- 
dition, since it scales time faster than space. When our results are applied to 
input processes satisfying a moderate deviation principle, the rate function 
of a multi-dimensional reflected Gaussian process appears. This establishes 
a rigorous link between large deviations of the behavior of heavily loaded 
queueing networks and those of reflected Gaussian processes. 

In this work we develop in detail this kind of relation for renewal arrival 
and service time processes. Our example partially parallels one of Puhal- 
skii [32] and yields the rate function of a reflected Brownian motion which 
similarly appears in the weak heavy traffic approximation of Peterson [29]. 
In a companion work [22] we apply the results of this work to a two queue 
network with long-range dependent traffic. There, large and moderate devia- 
tion principles and explicit logarithmic tail asymptotics for the steady-state 
queue length at the second queue are obtained. 

Our model captures the doubly infinite time interval. This permits us to 
treat directly steady-state situations. However, this feature requires the use 
of a topology which is slightly stronger than the usually considered topol- 
ogy of uniform convergence on compacts. In particular, the large deviation 
principle for the centered input processes must be provided in this topology 
before the main results of this work can be applied. By restricting the model 
to the positive time interval and by having initially empty queues, such a 
strengthening can be avoided. In addition, the considered model allows for 
mixtures of discrete customers and liquid inventory. 

Our construction of the network behavior on the doubly infinite time 
interval works without a definition of the state of the network, which is 
not obvious for multiclass fifo queues containing discrete customers and 
fluid. Hence, we cannot deal with networks starting from a fixed or random 
nonzero state at time 0. By modifying the input processes on the negative 
time interval, one can force the network into certain nonzero initial states, 
but a detailed treatment of nonzero initial states goes beyond the scope of 
this work. We therefore refrain from claiming uniqueness of the stationary 
network behavior and refer to [5] for results in this direction. 

An overview of this work is as follows: We recall necessary basic definitions 
and facts in Section 2. We state the main result for a single multiclass queue 
in Section 3. This result is extended to feedforward queueing networks in 
Section 4. We develop an application to renewal arrival and service time 
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processes in Section 5. Appendix A establishes moderate deviation principles 
for renewal processes in the topology considered here. Appendix B presents 
conditions for stationarity of network behavior. 

2. Preliminaries. A [0, oo]-valued lower semicontinuous function on a 
topological Hausdorff space is called rate function. A rate function is good 
if it has compact level sets. A sequence (Xfc)fcgN of random elements with 
values in a measurable space {E, £) satisfies a large deviation principle with 
normalizing sequence ihk)k&i and rate function / in the Hausdorff topology 
T on E \i limfc^oo bk = oo and for every measurable set Ag £, 

limsup-^logP(Xfc e A) <- inf I{x) 

fc-^oo Ofc ^^^^ 

and 

liminf -3-logP(Xfc G ^) > - inf I{x), 

fc— >oo Ofc x£A° 

where A^ (resp. A°) is the closure (resp. interior) of A in the topology T, and 
P is the underlying probability measure. For background on large deviations 
theory, we refer to [7]. See also [12] for background on large deviations of 
queueing systems. In this work the random elements are continuous time 
processes and the space E contains their sample paths which are introduced 
next. 

We let P be the set of functions d : M — > M which are right continuous, 
possess left-hand limits and have limits 



d:= lim d{t)/t 



and 



lim d{t)/t 



in M. On product function spaces these limits are taken componentwise. 
With I (resp. C), we denote the subset of P consisting of nondecreasing 
(resp. continuous) functions. We provide these (and further) function spaces 
with the cr-algebras generated by the family of one-dimensional projections, 
and topology induced by the norm || • || given by 

tm J- + \t\ 

With the usual pointwise addition and scalar multiplication, the spaces P 
and C become topological vector spaces. Product spaces are equipped with 
corresponding product cr-algebras and topologies. We shall mostly use bold 
face symbols for elements of, or mappings into, function spaces. 
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We let Pq be the subset of T> containing all functions d £T> satisfying 
d = d = 0. Clearly, if a function d G P satisfies d = d = : 5, then the centered 
function d — 5id is an element of "Dq, where id:M ^ R, 1 1— > t denotes the 
identity map of M. 

We recall that the composition o : (C x P) U ("D x X) — > P, (c, d) i-^ cod is 
continuous at elements in C x P and satisfies cod = cd and cod = cd (see 
Lemma B.l in [19]). 

On the set Pgup := {d G P:d > 0} we consider the "running supremum" 
sup : Psup Psup n I defined for d £ Pgup and t E M by 

(supd)(t):= sup d(r). 

Te]-oo,t] 

This mapping is continuous (see [11]), and satisfies supd = d and supd = 
max{0,d}. (For a finite subset IC of M, we let max/C denote the maximum 
of its elements.) 

On the set Xjnv := {c G XnC : c > 0,c > 0} we define the "right continuous 
inverse" Xinv ^X for c G Xinv and t G R by 

{c~^){t) := sup r. 

T£R,c{T)<t 

Clearly, c G Xinv implies coc"-*- = id < c'-'-oc, c~^ = 1/c and c~^ = 1/c. Here 
inequalities with vectors or functions are satisfied if they hold for every 
component or argument. 

For sets S and IC, we let 5^ be a set of (column) vectors with indices in 
/C and components in S. In particular, can be identified with S^^''"'''^ 
for A; G N. For v G and fc G /C, we let G 5 be the component of v with 
index k £ IC. 

3. Large deviations for a single multiclass queue. In this section we 
model the behavior of a queueing node with infinite buffer capacity and 
multiple customer classes. We state a large deviation principle for the be- 
havior of the queueing processes when the critical load 1 is approached, 
provided the input processes, which characterize arrival and service times, 
satisfy an appropriate sample path large deviation principle. 

We let be a finite set of customer classes. Two disjoint subsets 7i and 
C of the set M specify the high- and low-priority customer classes which are 
to be served by the node. Customers of a class in M.\{TCL) C) do not have 
to visit the considered node and leave immediately. Such classes are allowed 
in order to facilitate the extension to feedforward networks in the following 
Section 4. Customers of a class in 7i (resp. C) receive high-priority (resp. low- 
priority) service. Whenever a high-priority customer is present at the node, 
the service of low-priority customers is interrupted, and the customers with 
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high-priority are served using a fifo service disciphne, before service of low- 
priority customers is continued. Because of this preemptive-resume priority 
service disciphne, the high-priority customers experience a service by the 
queueing node as if there were no low-priority customers. Customers with a 
class in C share the remaining service capacity of the queueing node using 
a fifo service discipline. We assume that C is nonempty. This still permits 
modeling a pure fifo service discipline (without priorities) at the node by 
letting TC be empty. 

The external arrivals at the node are specified through elements a of 
I-^. The value of the jth component of such a nondecreasing arrival vector 
function a at time t is interpreted as the cumulative number of class j arrivals 
up to time t. This means that SLj{t) — aj{s) > counts the number of arrivals 
of class j customers during the time interval ]s,t] for s <t £M and j £ Ai. 
Here and in the following cumulative number of customers can be fractional 
(in order to permit scaling and to describe movements of fluid) and negative 
(in order to extend the cumulative customer numbers to the infinite negative 
past). Since the components of a are nondecreasing, there are no "negative" 
customers in the system. A jump in an arrival function aj can be interpreted 
as the arrival of a discrete customer with a size corresponding to the height 
of the jump, whereas a continuous increase corresponds to a movement of 
fluid. The vectors a,a G M.^ represent the asymptotical arrival rates on the 
negative and positive time intervals, respectively. 

The service time requirements of the customers needing service at the 
node are determined through a vector function s G I^!^^. For a customer 
class j GTiU C, the value Sj(6) specifies the cumulative number of served 
customers of class j when the node dedicates the total cumulative busy 
time 6 S M uniquely to this class of customers. In particular, the function 
sj^ (resp. s~''"oaj) describes the cumulative busy time required to serve the 
class j customers as a function of the cumulative number of arrived class j 
customers (resp. cumulative elapsed time). The continuity of Sj implies that 
the service time for a strictly positive number of class j customers is strictly 
positive. Jumps in the function sj^ correspond to a service interruption 
for class j customers. The vectors s,sgM^'-^'^ represent the asymptotical 
potential service rates on the negative and positive time intervals. 

In the following queueing model the vector function s only appears in 
terms of the form s~"'"oaj for j G A4. Hence, only the values of sJ"*" on the 
image {aj(t) :t € M} of aj are relevant for the behavior of the system. When 
the function aj ior j gTCU £. is piecewise constant, the sequence of jump 
sizes of the function sj^oaj can be interpreted as the service times required 
by the discrete customers modeled through aj. The behavior of the strictly 
increasing function s""*" between these jumps is irrelevant as long as the 
image of aj does not change. In Section 5 we will use these facts to model 
the vector functions a and s from sequences of inter-arrival and service times. 
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Pairs (a, s) G x satisfying the stability condition 

(2) E 

are called regular node primitives. The stability condition prevents overload 
during the infinite past, but does not forbid overload during finite time 
intervals or the infinite future. Given a regular node primitives pair (a, s), 
the cumulative idle time development of the node after service of only high- 
priority customers is defined as 

(3) U^(a, s) := sup ( id - ^ s-^oa^ | G J 

V jew / 



The (backlogged) workload behavior of high-priority customers (i.e., the 
time needed to complete service of all high-priority customers which are 
currently at the node) can be expressed as 

(4) V^(a, s) := ^ s-^oaj - id + U^(a,s) G V. 

For Ti = (pure fifo service discipline), we get U^(a, s) = id and V^(a, s) = 
0. 

The remaining unused potential busy time of the node is dedicated to 
serving low-priority customers. Therefore, the cumulative idle time develop- 
ment of the node is defined by 



Y^'^(a,s) :=sup( U^(a,s) - E^7^°^i ) ^ ^• 



(5) 

The (backlogged) workload behavior of low-priority customers is given by 
(6) W^'^(a,s) := ^ sj^oaj - U^(a, s) + Y^'^(a, s) G V. 

Lemma 2 of [21] shows that these definitions imply 

V^(a,s)>0, 
V''(a,s) = 0, 

/ V^(a,s)(i)dU^(a,s)(t) = 0, 

W^(a,s) > 0, 
W^(a,s) = 0, 



/ W^'^(a,s)(t)(iY^'^(a,s)(t) = 0, 

Jm. 

f V^(a,s)(t) dY'^'^{a,s){t) = 0. 
Jr 
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Here du denotes the cr-finite measure induced by a nondecreasing path u G T 
on M. Conversely, this lemma also states that these properties in conjunction 
with equations (2), (4) and (6) imply equations (3) and (5). Hence, the 
definitions in equations (3) to (6) are justified by the nonidling service, the 
preemptive-resume priority structure of the service discipline, the wish for 
stability on the negative time interval, and the interpretation of s""*" o as 
the cumulative busy time required to serve all arrivals of class j customers 
as a function of time. 

Lemma 2 of [21] also establishes the hidden identities 



Y^'^(a,s) =sup( id- ^ sT^oa 



jenuc / 



W^'^(a,s) + V^(a,s) = ^ sj^oaj - id + Y^'^(a, s). 

jeHuc 

Next, we define the cumulative departures D^''^(a, s) gX-^ at the node. 
The jth component in this vector of functions describes the cumulative 
number of customers of class j £ M which depart from the node as a function 
of time. We define it for t G M by 



sup (■'")) if i £ 



sup (''")) if j £ -C, 

^aj(t), otherwise. 

The conditions below the supremum in this definition assures that the work- 
load incorporated in high- or low-priority customers at the node reaches at 
least its current workload. Due to partially served customers, it can exceed 
the current workload. Customers of each priority group leave in the order of 
their arrivals, and all customers of one priority group arriving at the same 
time instant leave together. Conversely, a larger or smaller number of de- 
partures would not be consistent with this order of the departures within 
the priority groups, the requirement that served customers must leave, or 
the definition of the workload behavior. In this way Y)"^'^ models fifo service 
discipline within the class of high- or low-priority customers. Customers of 
classes which need not visit the node leave immediately. 

The behavior of queue lengths can now be defined as the difference be- 
tween cumulative arrivals and departures 



Q^'^(a, s) := a - D^'^(a, s) G V 



M 
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Finally, the time at which an observer (e.g., an additional customer which 
does not affect the behavior of the system) of class j £ M entering the node 
at time t €M leaves the node is 

{inf T, ifjen, 

T>t, D"'^(r)>af(t) for every ^GW 
inf T, ifjGA 

T>t, D^'''{T)>at(t) for every ieC 
t, otherwise. 

'hi • A// 

Hence, Z^- ' (a, s) — id G represents the sojourn time of a class j observer 
at the node as a function of its arrival time. 

The model could be extended to incorporate more than two priorities; 
compare [4]. Since high-priority customers vanish in the heavy traffic limit 
(see Remark 4 after Theorem 1 below), we refrain from such an extension. 

Similar pathwise models are frequently used to analyze the behavior of 
queueing systems [4, 13]. The one presented here is special in that it unifies 
discrete customer and fluid models. This is particularly useful in our context 
since sample path large deviation scalings often lead to a focus on continuous 
functions (i.e., the rate function is infinite for noncontinuous arguments) 
corresponding to movements of fluid through the system, even when the 
cumulative arrival (and, hence, also the departure) processes model discrete 
customers (i.e., they are concentrated on piecewise constant functions). 

Our model includes the infinite negative time interval which is particularly 
advantageous in analyzing stationary behavior [11]. We refer to Lemma B.2 
in Appendix B for independence and stationary increments conditions which 
are sufficient for a stationary behavior of the queue. 

On the other hand, by restricting the elements of all considered function 
spaces to the time interval M+, and, furthermore, to those starting in 0, the 
model captures the behavior of the multiclass node when it starts empty at 
time 0. Clearly, in this modification, the supremum in the definitions of the 
mappings || • ||, sup, and has to be restricted to nonnegative arguments, 
and condition (2) becomes meaningless. 

The full model with time interval M and its restriction starting empty at 
time zero are intimately related: If (a, s) is a regular node primitives pair 
and i G M a time where the associated queue is empty [i.e., 'W'^'^{a.,s){t) = 
Y'^{a,s){t) = 0], then, for every t' >t, the supremum in the definitions of 
U'^(a,s)(t'), Y'^''^(a,s)(t') and D^"^(a, s)(t') can be restricted to the fi- 
nite interval [t,t']. This observation implies that if (a,s) is a random node 
primitives pair satisfying (a(0), s(0)) = and the assumptions of Lemma B.2 
(implying the stationarity of the associated behavior of the queue), work- 
load and queue length distributions in the model restricted to the positive 
time interval, converge to the corresponding stationary distributions as time 
increases to infinity. This property carries over to the network case of Sec- 
tion 4. 
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Next, we specify a time-homogeneous critical load situation: For the re- 
mainder of this section, we fix vectors a G M^^^'^ and a G M^* satisfying 



o->0, 



(7) /:=^^>0 
and 

If we interpret a as vector of arrival rates and a as vector of service rates, 
the last condition means that the node is critically loaded. Condition (7) is 
equivalent to the existence of a j € £ with aj > 0. 

In order to formulate the main theorem of this section, we define for 
functions a G T>-^ and s G j)HyJC satisfying 



jeHuc 



the functions 



and 

jeHuc 

Furthermore, we let G be the vector with components 



^^'^(a,s):= y ^-^opiiid +Y«'^(a,s)GP. 



a 



[O, ifjGAl\£, 
and e"^ the vector in with ej' = 1, if j G >C, and ej' = 0, otherwise. 

Theorem 1. We let {aLk)km be a sequence of processes on ana 
(sfc)A;eN CL sequence of processes on 1^1^^ . We assume that, for every A; G N, 
there are vectors a^. G M:^ and G W^'^^ such that, with probability one, 

ak=kk = a/c, 
o-fe = Sfc = Sfc > 0, 

(10) Pk:= E ^<1- 

jeHuc 
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We assume that the sequence of centered regular node primitives processes 
(a/c - afcid,Sfc - o-fcid)^ 



satisfies a sample path large deviation principle with normalizing sequence 
(&fc)fceN o,nd good rate function I on the space assume 
that I takes the value oo whenever a component of its argument functions 
is not continuous. 

If there is a sequence {dk)keN 'i'n IR+ and vectors a E R-'^, a G M^^'-- with 

(11) oo = lim dk, 

(12) a = lim (a^, — d^ctt), 

fc— >oo 

(13) a = lim {a^ - dkcr) 

K— >00 

and 

(14) ~p:= E (^^-^)>0, 

the following 8 statements are valid with normalizing sequence {hk)k&n- 
1. The sequence {dk\J'^{ak,Sk) - dk{l - p^)id)fcgN, where 



Pk 



satisfies a large deviation principle on Dq with good rate function 
ui— > inf /(a, s), 

a£V^,seT>y-^'^, u=U'W(a+aid,s+5id)-p'Wid 

where 

pT^ ■=Y 



2. The sequence (dfcV^(afc, Sfc))fcgN satisfies a large deviation principle on 
T>Q with good rate function 

0, i/v = 0, 
oo, otherwise. 

3. The sequence (dfcY^"^(afc, s^) — dk{1 — Pk)'i^)k& satisfies a large devia- 
tion principle on Dq with good rate function 

y I— > inf /(a,s). 

aGX'^,se-D5^u£^y^Y«'^(a+aid,s+5-id)-pid 



12 K. MAJEWSKI 

4. The sequence ((ifcW^''^(a,t) Sfc))yfcgN satisfies a large deviation principle on 
Do with good rate function 

w h->- inf /(a, s). 

aeO(t^,sGX'2^^^,w=W«.^(a+aid,s+5id) 

5. The sequence (D^''^(afc, s^) — afcid)/;^^ satisfies a large deviation princi- 
ple on Vq^ with good rate function 

d I— > inf /(a, s). 

a£V^ .sGXlJ^ur ,d=a-a£ W'W-C (a+aid,s+5-id) 

6. The sequence (Q^''^(afc, Sfc))fcgi^ satisfies a large deviation principle on 
T>Q^ with good rate function 

qH-> inf_ /(a, s). 

aeV^ .sef^^^^ ,q=a^ W«.-C (a+aid,s+5-id) 

7. The sequence ((ifcZ^'-^(afc, s^) — (ifcid)fcgN satisfies a large deviation prin- 
ciple on T)q^ with good rate function 

z 1-^ inf /(a, s). 

aeV^ ,sGX55^^^ ,z=e': W«'£ (a+aid,s+5-id) 

8. T/ie previous statements are also valid in combination: For instance, 
the sequence ((ifcW^''^(afc, s/;), Q^'^(afc, Sfc))fcgpj satisfies a large devia- 
tion principle on Vq x Dq^ with good rate function 

(w,q)H^ inf I(a, s). 

a e D,-^ , s e Vy-^'^,V! = W'«'-C(a + aid, s + aid), 
q = a^'W^'^(aL + aid, s + aid) 

The statements of the theorem remain valid when the elements in all un- 
derlying function spaces are restricted to the time interval M_|„ and to those 
starting in 0, and the definitions of the model are modified accordingly. 
[Hence, conditions (2) and (9) become meaningless] With this modification, 
the statements are also true when conditions (10) and (14) are dropped, 
and/ or the topology is changed to the topology of uniform convergence on 
compacts. 

Before giving the proof we remark: 

Remark 1. Limits (11) to (13) imply 

lim — = a, 

fc— >oo dk 

lim — = o", 

k-^oo dfc 

lim pk = l 

k^oo 
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and 

hm 4(1 -pfc) = hm > — i J-^ = p. 

Hence, the asymptotic load at the queue on large time intervals must ap- 
proach the critical value 1 with a speed of order 1/4 as A; — > oo. This prereq- 
uisite is the reason for speaking about large deviations "in critical loading" 
as in [32]. 

Remark 2. The requirement (14) is not as restrictive as it might appear 
since conditions (10) to (13) imply /5 > 0. 

Remark 3. The mappings U^, Y^'"^ and W'^"^ appearing in the ex- 
pressions of the rate functions for the network behavior have a much simpler 
structure than the network mappings (particularly the mapping Yf^''-' incor- 
porating the queueing discipline). They consist mainly of linear combinations 
of the argument functions and the sup-mapping. This reflects the fact that 
the behavior of the queueing node simplifies close to a time-homogeneous 
critical load; compare [21]. 

Remark 4. In statement 8 of the theorem the rate function for workload 
and queue length behavior is only finite when the ith component of the queue 
length behavior is af times the workload behavior. In this case, high-priority 
customers do not experience a delay. An analogous state space collapse is 
well known from heavy traffic convergence in distribution [29]. 

Proof. We first prove statement 4. Under the assumptions of the the- 
orem, we define for every E N the mapping : x V]^^^ Vq by 

F.(a,s):=4W^.-(^±^,i±^ 
V dk dk 

if a e and s G P^^^ are such that a+ a^id G J-^ and s + a^id e I^';'^. 
Otherwise, we set 

Ffc(a, s) := W'^'^(a + aid, s + crid). 

Conditions (10) and (14) ensure that the mapping F^ is well defined. 

We let a G T>^ and s £ V'^^^ with /(a, s) < cx3. Hence, a G C-^ and s G 
qHvjc bgcausg J(a, s) would be infinite if a component of a or s were not 
continuous. If (afc)^^^ is a sequence of functions in Vq^ converging to a and 
(sfe)fcgN is a sequence of functions in D^-^^ converging to s, then 

(15) lim Ffc(afc, Sfc) = W''^"^(a + aid, s + aid). 
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This is true if + Qfcid G and + cifcid € ZT^^ for every /c G N because 
of the fohowing: In this case we set, for /c G N, 

afc + Offcid 
afe := ^ 

and 

Sfc + (Tfcid 
Sfc := 1 , 

dk 

and thus obtain a sequence of regular node primitives (afc,Sfc)fcgp^ which 
satisfies 

hm (ifc(afc — aid) = hm (a^ + (a^ — (ifca)id) = a + aid 

fe— >00 fc— KX> 

and 

hm dk{sk - o"id) = lim (sfc + (dfc - (ifccr)id) = s + aid, 

fe— >oo fc^oo 

because of conditions (12) and (13). Thus, this sequence satisfies the pre- 
requisites of Theorem 4 in [21] and we obtain 

hm Ffc(afc,Sfc)= hm 4W^'^(afc, s^) = W^'^(a + aid, s + aid). 

fe— »00 K— >00 

If afc + afc ^ or Sfc + ^ ^^fv ^ for some G N, we stiU obtain the 
convergence (15) since W^''^ is continuous. 

We can therefore apply the version of the contraction principle of [31] (see 
Remarks 2 and 3 after Theorem 2.1 therein) in order to conclude that the 
sequence Ffc(afc — Ofcid, Sfc — tifcid) satisfies a sample path large deviation 
principle with the normalizing sequence {bk)keN and the good rate function 
given in statement 4. Statement 4 now follows from the fact that, for every 
A; G N, we have, with probability one, 

Ffc(afc - ttfcid, Sfc - o-fcid) = dfcW^'^ ( £^ ) = 4W^'^(afc, Sfc). 

\dk dkJ 

In the last equation we used the property of the mapping W^'^ that for 
every regular node primitives tuple (a, s) and value ^ > 0, we have 

W^'^ (^a, ^s) = W'^'^ (a, s) . 

The proofs of the other large deviation principles are analogous. However, 
in certain cases the previous invariance must be replaced by the positive 
homogeneity as exemplified by the equation 

D«'^(ea,es)=eD^'^(a,s). 

This distinction is responsible for the fact that in statements 5, 6 and 8 the 
factor dfc disappears. 
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When the model is restricted to the time interval and to functions 
starting at 0, the proofs are completely analog and do not make use of 
conditions (10) and (14). Compare Remarks 5 and 6 after Theorem 1 in [21]. 
In this situation the behavior of the queue on the time interval [0, t] for t > 
depends only on the cumulative arrivals until time t and the cumulative 
number of served customers up to the cumulative busy time t. This permits 
us to transfer the statements to the topology of uniform convergence on 
compacts and completes the proof. □ 

4. Large deviations for feedforward networks. In this section we consider 
n G N queues numbered from 1 to n with fifo and priority service disciplines. 
Therefore, several symbols from the previous section receive an additional 
index i specifying the queue with which they are associated. Each customer 
of the network belongs to a class of the finite set of customer classes Ai. 
For every queue i G {l,...,n}, disjoint subsets 7ii and Ci of Ai specify 
the classes of customers which must visit the zth queue. Before leaving the 
network, the customers of class j € 7W have to visit all queues i < n with 
j GfiiU Ci in increasing order. Because of this last requirement, our model 
only captures feedforward networks. 

Customers of a class in 7ii receive high-priority and customers of a class 
in Ci low-priority service at queue i < n. In particular, the priority group to 
which a customer class belongs can change from station to station. When 
high-priority customers are present at a queue, the service of low-priority 
customers is interrupted and then resumed when no high-priority customers 
require service. Among the customers of one priority group the service time 
is shared in fifo order. We assume that Ci is nonempty for every i < n. 
By having Tii = for some or all nodes i < n, one can model stations or 
networks with a pure fifo service discipline. 

As for the single node, external arrivals at the network are specified 
through elements a of T-^ . The service time requirements of the customers 
at the queues they have to visit are determined through vectors of functions 
Si G Xjj^^ \ . . . , s„ G T;^" " . For every queue i <n and customer class 
j (z TCiU Ci, the value of the component Sjj at b specifies the cumulative 
number of class j customers which is served by queue i when it dedicates 
the cumulative busy time b to this class of customers. 

In order to be able to construct the behavior of the queueing network for 
given cumulative arrival function vector a and service time function vectors 
si , . . . , s„ , we have to guarantee that the queues have not been overloaded 
during the infinite past. This amounts to requiring the stability condition 



(16) 




a, 



for i = 1, . . . , 



n. 
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An (n + l)-tuple of function vectors (a, Si , . . . , s„) with a G , Si G I^^'^^'^ , 
. . . , s„ E ^Tj^"^"^" satisfying this condition is cahed a regular network prim- 
itives tuple. For such a regular network primitives tuple (a, Si, . . . , s„), we 
define 

Do(a,si, . . . ,s.„) := a 
and, inductively, for i = 1, . . . , n, 



U,(a, 


Si,. 


• ) ^n) 


:=U^'(Di_i 


a,si, . . . 




), 


Vi(a, 


Si,. 


• ) Sn) 


:=V^'(Di_i 


a,si, . . . 




), 


Yi(a, 


Si,. 


• ) Sji) 


:= Y^''^»(D^ 


_i(a,si. 


• • ) ) 




Wi(a, 


Si,. 


• ) Sn) 


:= W'^"^»(D, 


_i(a,si, 


• • • ) Sriy 


,s 


D,(a, 


Si,. 


• ) Sji) 




_i(a,si, . 


• • 5 ) 




Qi(a, 


Sl,. 


• ) Sn) 




_i(a,si, . 


• • J ) 




Zi(a 


Si, . 


• ) Sn) 


:= Z^'''^'(Di_ 


_i(a,si, . 


• • ) Sn) ) 





and 



These mappings specify the behavior of the zth node in the network anal- 
ogously as in the single queue case. Again the model contains stationary 
(see Appendix B) and initially empty queueing networks (by restricting the 
functions in the underlying function spaces to the time interval 1R+ and to 
those starting in 0, and modifying the definitions of the model accordingly) 
as special cases. In the latter case condition (16) becomes meaningless. 

satisfy, for every i<n, 



We let vectors a£R^, ai G R^^^^\ . . . , (t„ G M^„u£„ g-^^^ ^j^-^j^ 



(17) 

and 
(18) 



o-,; > 0, 



E — = 1 



When Oj is interpreted as the exogenous arrival rate of class j GTCiU Ci 
customers and (Tjj- as the service rate of class j customers at queue i <n, 
the last condition says that node i is critically loaded. 

We shall need the vector af G M-'^ for i <n, which is defined by 



(19) 



ifiGA^\£.. 
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Furthermore, we define the strictly lower triangular matrix G G M"^" by 
(20) Gi,h : 



^ Jhl^ ifh<i, 



. 0, otherwise. 

The triangularity of the matrix G is a consequence of the feedforward struc- 
ture of the modeled network. 

Close to the "critical" linear network primitives tuple (aid, diid, . . . , cr„id), 
the behavior of the network simplifies [21]. In order to model the simplified 
behavior we introduce for i = 1, . . . ,n the P- valued mappings Xj, Wj, Uj 
and the X- valued mapping Yj by setting 



X.(a,s„...,s„):= E (^-fl^Iof^)), 



Yi{a, si, . . . , s„) := sup ^-Xi(a, si, . . . , s„) + ^ Gi,/,Wi(a, si, . . . , s„)^ , 
Wi(a, si, . . . , s„) := Xj(a, si, . . . , s„) 

i-l 

)+Yi(a 

h=l 

U. (a. s.. .... s„) := i: f (gg) - ^ + g "fe^^'-'"- . 

These mappings are well defined for argument function vectors a G T)^ , 
si G P^iU^i, . . . , s„ G pWnU£„ satisfying 

(21) J2 (=^-^-^)>0 for every i<n. 

Theorem 2. We let (a/c)fcgf^ be a sequence of cumulative arrival pro- 
cesses in and {sk^i)k(m o, sequence of cumulative service time processes 
inX'^'-'^^' for every i <n. We assume that, for every A; G N, there are vectors 
oik G CFk,i G . . . , Ok,n G M^"^^" such that, almost surely, 

ak=tk = a/c, 

= ifc,i = Sfc,i > fori = l,...,n, 

(22) pk,i:= fori = l,...,n. 
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We assume that the sequence of centered regular network primitives tuples 

(afc - a/fcid,Sfc^i - (TA;,iid, . . . ,Sfe^„ - afc,nid)fc6N 

satisfies a large deviation principle with normalizing sequence {bk)keN o-i^d 
rate function I on x V^^^'^^ x ••• x assume that 



this rate function takes the value oo whenever a component of its argument 
functions is not continuous. 

If there is a sequence (4)fceN in ond vectors a G IR-^, ai G ro^i'-'-^i 





with 




OO 


= lim dk, 

fe— >oo 




Q 


= lim (dfc 

fe— >oo 


- dko), 




= lim ((Tfc,, 

fe— »oo 


; - dkUi) 


Pi ■■ 


- s ( 

jeHiUCi ^ 


aij_aj_ 



fori = l,...,n, 

(23) pr.= i: (?^3-|^)>0 /or.= l,...,n, 

the following 8 statements are valid with normalizing sequence {bk)ken' 

1. For every i < n, the sequence {dkUi{ak,Sk,i, ■ ■ ■ ,h,n)-dk{l- Pk^i)id)k(zn, 
where 



satisfies a large deviation principle on Dq with good rate function 

uh^ inf /(a,si,...,s„), 

& G , SI e , . . . , s„ e T^J'-u^" , 

u = Ui(a + aid,si + aiid, . . . , s„ + a-„id) — p^id 



where 



2. For every i <n, the sequence (dfcVj(afc, i, . . . , Sk^n))k&i satisfies a large 
deviation principle on Dq with good rate function 

0, i/v = 0, 
oo, otherwise. 

3. For every i < n, the sequence (4Yj(afc, Sfc,i, . . .,Sk,n) - 4(1 - Pk)id)k 
satisfies a large deviation principle on Dq with good rate function 



inf /(a, si 

y = Yi(a + aid,si + Slid, . . . ,s„ + o-„id) - ftid 



) ■ • ■ ) ^nj 
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4. For every i<n, the sequence ((ifcWj(afc, s^^i, . . . , Sk^n))k&^ satisfies a large 
deviation principle on Dq with good rate function 

w 1-^ inf /(a, §1, . . . , s„). 

w = Wi(a + aid, si + aiid, . . . , s„ + CTnid) 

5. For every i <n, the sequence (Dj(a,fc, s^^i, . . . , Sfc^„) — afcid)fcgN satisfies a 
large deviation principle on Dq^ with good rate function 

inf 7(a,si, . . . ,s„). 

d = a — a^Wi(a + aid, si + CTiid, . . . , s„ + (Tnid) 

6. For every i<n, the sequence {Q,i{sLk,Sk^i, ■ ■ ■ ,Sk^n))k€N satisfies a large 
deviation principle on Vq^ with good rate function 



qi-^ inf /(a,si,...,s„). 

q = af Wi(a + aid,si + aiid, . . . , s„ + a-„id) 



7. For every i <n, the sequence ((ifcZj(afc, i, . . . , „) — (ifcid)fcgp^ satisfies 
a large deviation principle on "Dq^ with good rate function 

inf /(a,si,...,s„). 



,HiU£i 

e-''Wi(a + aid,si + (Tiid, . . . ,s,i + 5-,iid) 



8. The previous statements are also valid in combination: For instance, for 
every i,h<n, the sequence 

(dfcWj(afc,Sfc^i, . . . ,Sk^n), Q/i(afc,Sfc,i, . . ■,Sk,n))ken 

satisfies a large deviation principle on Vq x Pg^ with good rate function 

(w,q)H^ inf /(a,si, . . . ,s„). 

w = Wi(a + aid, si + o"iid, . . . , s„ + 5-„id), 
q = a^W|j(a + tiid, si + 5"iid, . . . , s„ + (T„id) 

The statements of the theorem remain valid when the elements of all un- 
derlying function spaces are restricted to the time interval and to those 
starting in 0, and the definitions of the model are modified accordingly. 
[Hence, conditions (16) and (21) become meaningless.] With this modifica- 
tion, the statements are also true when conditions (22) and (23) are dropped, 
and/or the topology is changed to the topology of uniform convergence on 
compacts. 



The remarks on Theorem 1 similarly apply to each queue individually. 
We make the following additional remarks concerning this theorem: 
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Remark 5. We say that a pair (w,y) G 2?" x solves the Skorokhod 
problem for the vector z E D" and the matrix R G M"^" (on the time interval 
M) when the following conditions are met for every component z < n: Wj > 0, 
Wj = 0, Wi = Zj + X;h:=i Ri,hyh, and /jj vifi(t) dyi(t) = 0; compare [9, 14, 20]. 

Here we define the matrix R as the inverse of the matrix I + G, where / is 
the n X n-identity matrix. Due to the feedforward structure of the queueing 
networks considered in this work, the matrix R is triangular and the matrix 
consisting of the absolute values of the entries of the matrix I — R has 
spectral radius 0. Hence, Theorem 6 in [20] states that, for every z in the 
set 

7^:={zGP": (/ + G)z < 0} 

(with the usual definition of a matrix- vector product), the Skorokhod prob- 
lem for z and R has a unique solution which depends continuously on zGTZ. 
We denote the map from z G 7^ to the vector of functions w G in the 
unique solution (w,y) of the Skorokhod problem for z and R with <I>. Such 
a solution map is called a Skorokhod or reflection map. 

By Lemma 1 of [21], we have W(a, si, . . . ,s„ ) = ^>(i?X(a,si,...,s„)) for 
all a G V-^, si G P^i^^i^ . . . , s.„ G V'^"^'^^ satisfying condition (22). Under 
the assumptions of the theorem, we therefore obtain that the sequence of 
n-dimensional processes ((ifcW(afc, s^^i, . . . .,Sk,n))k&i satisfies a sample path 
large deviation principle in Pq with good rate function 

inf /(a,si,...,s„). 

w = ^{RX-ia + oid, si + triid, ...,§„ + 5-nid)) 

In this light the rate function for the sequence of n-dimensional workload 
processes in the theorem can be looked upon as a "reflected" rate function in 
the sense that it is derived from the original rate function through a reflection 
map. Similarly, one can represent the rate function for queue length and 
sojourn time processes as "reflected" rate functions. 

Remark 6. When the rate function / in this theorem stems from a 
moderate deviation principle, it is usually also the rate function in a large 
deviation principle for the tail probabilities of a Gaussian process. In view 
of the previous remark and of the fact that the reflection map is positively 
homogeneous [i.e., <I>(cz) = c<I>(z) for every c > and z G 7^], the rate func- 
tion for the sequence of n-dimensional workload processes in this case is 
also the rate function in a large deviation principle for the tail probabili- 
ties of an n-dimensional reflected Gaussian process. Similar statements hold 
for the queue length and sojourn time processes. In this sense our work 
can rigorously justify the analysis of large deviations of reflected Gaussian 
processes [1, 20, 23, 25, 27, 28] as a means to analyzing logarithmic tail 
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asymptotics of the behavior of queueing systems in critical loading. Com- 
pare Section 5 for an example with renewal and [22] for an example with 
long-range dependent input processes. 

Remark 7. This theorem parallels the convergence of the queueing pro- 
cesses in distribution to reflected (fractional) Brownian motion in heavy 
traffic, when the sequence of random network primitives tuples satisfies a 
functional (fractional) central limit theorem [21, 29]. 

Remark 8. Statement 7 of Theorem 2 reveals the validity of the "snap- 
shot principle" [34] in large deviation asymptotics. Namely, the rate function 
for the sojourn time behavior is only finite, when at every time point the so- 
journ time of an observer arriving at this time is a specific linear combination 
of the workload of the nodes at this time. 



Proof of Theorem 2. This is analogous to the proof of Theorem 1, 
but using Theorem 5 of [21] instead of Theorem 4 of [21]. Due to the feedfor- 
ward structure of the network, one can also prove this theorem by iteratively 
applying Theorem 1. □ 

5. Example: renewal input processes. In this section we show how The- 
orem 2 can be applied in a case where the inter arrival and service times 
of the customers of each class form independent and identically distributed 
sequences of random variables. We partially parallel Theorem 4.1(c) of [32], 
but since we aim for the stationary case, and the theory of sample path large 
deviations in the topology || • || is less developed, we need stronger condi- 
tions. We show that the obtained rate function is in fact the rate function in 
a large deviation principle for scaled reflected Brownian motion. When at- 
tention is restricted to the positive time interval and initially empty queues, 
one can use analogous conditions as in Theorem 4.1(c) of [32] and mimic the 
development of this section in order to obtain large deviation principles for 
the queueing behavior in the topology of uniform convergence on compacts. 

Given a sequence {Xfi)hGZ of strictly positive, identically distributed and 
independent random variables with finite mean, we define the linearly in- 
terpolated partial sums (lips) process X^'p^ = (X^'P^(t))tgK by setting, for 

0, if/i = 0, 

h 



(24) Xi^P^(/i) := 



J2 Xra, if /l > 0, 

=1 

if/i<0, 



m=l 



m=l 
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and linearly interpolating these values between subsequent integer time 
point s. Through the law of large numbers we get X'^P*^ G C n J and X^'P" = 
X'^P*^ = E(Xi), with probability one. By modifying the underlying probabil- 
ity space on a zero-set, we can and will assume without loss of generality 
that these properties are satisfied for all sample paths of linearly interpolated 
partial sums processes in this work. 

Clearly, the time shifted process X^'p^(- + h) — X''P'^(/i) has the same dis- 
tribution as X'^P*^ for every /i G Z. Furthermore, the process (X^'p**)"-*- is a 
renewal process having a jump at time and being linearly interpolated 
between subsequent jump times. 

For X G M, we let [x\ be the largest integer h £Z with h < x. li N is an 
uniformly distributed random variable on [0, 1[ independent of (-^h)/iGZ; the 
process x on T defined by x(t) := [(X^'p** — NXi)~'^{t)\ for t G M is a renewal 
process with stationary increments [i.e., the distribution of x(- +t) — x(t) is 
the same as the one of x for every t G M]; see Exercise 1.2.2 in [3]. 

In the following we shall use these observations to define arrival and ser- 
vice time processes from sequences of random variables. We scale time of the 
kth input process by the factor k and its state by l/\/bkk, where {bk)keN is 
a sequence of positive numbers satisfying 

(25) lim ^ = 
and 

(26) lim = oo. 

fc-»oo log k 

We first specify the sequence of renewal arrival processes: For every G N 
and J G Al, we let {Akj^h)hez be a sequence of identically distributed strictly 
positive random variables having finite second moments. For every A: G N and 
j ^ Ai, we let Nf^j be a random variable which is uniformly distributed on 
[0,1[. We assume that for every fixed A; G N the random variables in the 
set {Akj^h : i G A^, /i G Z} U {N^j : j G A^} are independent. We define the 
process on I-^ by setting, for j £ A4 and t G M, 

a,,,(t) := ^L(a5 - N,,A,,,,r\kt)\. 

Hence, the process a^j specifies discrete arrivals of class j customers with 
size 1/yJhkk having independent inter arrival times distributed as Ai^j^i/k. 
We define, for /c G N and j £ M, 




For every /c G N, the process has stationary increments and satisfies 
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Second, we specify service time processes: For every A; G N, i < n and 
j £7iiUCi, we let {Sk,ij,h)h<=z be a sequence of identically distributed strictly 
positive random variables having finite second moments. We assume that, 
for every fixed k gN, the random variables in the set {Sk^ij^h'-i 1^n,j G 
HiU £i, heZ}U {Akj^h -.j e M,heZ}U {Nkj -. j e M} are 'independent. 
We define the process s^^i on Tj^^^ ' by setting, for j GTiiU Ci and t G M, 

1 



The service time process j j specifies independent service times distributed 
as the random variable Sk,i,j,i/k for discrete class j G TLi U Ci customers of 
size l/\/bkk at queue i. We note that the linear interpolation involved in the 
definition of Sk^ij does not disrupt this interpretation, since in the definitions 
of the behavior of the queueing network the process s^^ij does not appear 
outside terms of the form s^]"^oafc j. 

We define, for /c G N, i < n 'and j GTiiUCi, 



1 



bk E{Sk,ij,i) 

For every /c G N, and i <n, the process Sk^i satisfies 

(^k,i — ^k,i — ^k,i' 

The load generated by the sequence of input processes must approach a 
critical load as k increases. This is captured by assuming the conditions 

(27) lim EiAk,j,i) = aj^ forjGTW, 

(28) lim E{Sk,ij,i) = cr^f for i < n and j GHiU Ci, 

where a G M^, ai G M+''^^\ . . . , cr„, G M^"U£„ "critical rates" vectors 
with strictly positive components and satisfy conditions (17) and (18). In 
order to reach the critical load at a specific speed, we additionally assume 
that there are vectors a G M-^, ai G M'^i^^^, . . . , (T„ G M'^"^^" satisfying 

for j G M, 




- (Ji.jj for i < n and j G HiU Ci. 

Next, we define a rate function for the sequence of input processes. We 
assume that there exist vectors u G M^, vi G R^^^^^ , . . . , t;„ G R^"^-^" with 

^lim Var(^fcj i) = ii| for j £ ^A, 

lim Var(5fc i j i) = vf a for i<n and j G U Ci. 
fe— »oo ' ' ' " 
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We define the good rate function /Brown : ^ j {cx)} for x G Pq by 

jBrown^^^ ■= ) — ^ — dt, if x(0) = 0, and x is absolutely continuous, 

[ oo, otherwise. 

Here x denotes a derivative of an absolutely continuous function x G T>q. 
This rate function is well known from Schilder's theorem [8] dealing with 
logarithmic tail asymptotics of Brownian motion. We let P'^'^^"^'^^ be the good 
rate function on x p^iU^li x . . . x pJ^-uA4„ ^^^^^ 

^SjBrown^;^^.) n ^3 jBrown^g.^^.) 



/~'(a,si,...,M:=E ' ,2 ' +E E — 



2 

II ■ ??- ■ 

Here we set x/0 := 0, if x = 0, and x/0 := cxo, if x G M \ {0}, in order to 
capture the case where some of the limiting variances are zero. We can now 
apply Corollaries A. 2 and A. 3 of Appendix A in order to get the following: 

Corollary 3. // there exists 5 > and c > such that, for every A; G N 
and y G [—5, 5], 

log£;(exp(2/(^fcj- 1 - E{Akj,i)))) < cy^ for j G M, 

logE{exp{y{Sk,i,jA - E{Sk,ij,i)))) < cy^ fori<n and j GTiiU Ci, 

the sequence 

(afc — afcid, s^^i - (7fc,iid, . . . , Sfc^„ — (Tfc,„id)fcgN 

satisfies a sample path large deviation principle on Vq^ x V^^'^'^^ x • • • x 
jyHnVJCn yjj^f^f^ Qood ratc function p^^™^^ and normalizing sequence 

If for a /c G N the stability condition (22) is satisfied, the random network 
primitives tuple (a^, i, . . . , Sfc^„) is regular, and the queueing processes 
Ui(afc,Sfc^i, . . . ,Sfc^„), Vj(afc,Sfc^i, . . . ,Sfc^„), Yj(afc, s^^i, . . . , Sfc^„), Wj(afc,Sfc^i, 
■ • • ,^k,n), Dj(afc,Sfc^i, . . . ,Sfc^„), Qi(afc,Sfc^i, . . . ,Sfc^„) and Zj(afc,Sfc^i, . . . ,Sfc^„) 
are well defined for every queue i <n. From Lemma B.2 of Appendix B, we 
deduce that the processes V(afc, s^^i, . . . , Sfc^„), W(afc, s^^i, . . . , Sfc^„), Q(afc, Sfc^i, 
. . . , Sk^n) and Z(afc, Sfc^i, . . . , Sk,n) - id are jointly stationary. 

Corollary 4. //, for every k £N, conditions (22) and (23), and the 
assumptions of Corollary 3 are satisfied, the prerequisites and consequences 
of Theorem 2 are valid with I : = p^^'^'^^^^ and d^ ■= \/kJb^ for /c G N. 

Clearly, independent exponentially or deterministically distributed inter- 
arrival and service times with mean values satisfying conditions (22), (23), 
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(27), (28), (29) and (30) constitute a possible choice through which the 
assumptions of the previous two corollaries are satisfied. 

For a positive semidefinite covariance matrix V G M*^^" and a vector b € 
we define the good rate function ly ■ T>q ^ R U {oo} by 



^ dt, if X is absolutely continuous 

and x(0) = 0, 
otherwise. 



Here x(t) is the vector of derivatives of the absolutely continuous function 
X € at time t G M, and x(t)^ its transpose. If the matrix V is degener- 
ate, the expression x(t)-'"V~-^x(t)/2 in this definition must be understood 
as supj^gign(y"^x(t) — y^Vy /2)\ compare Remark 4.1 in [32]. We note that if 
B is an n-dimensional Brownian motion with stationary increments on the 
time interval M, passing through at time 0, having covariance matrix V 
at time 1, and possessing drift 0, the sequence ij^ / y/bk)k<^^ satisfies a large 
deviation principle with normalizing sequence {bkjk^-^ good rate function 

Here we define V := RUR^ , where is the transpose of the matrix R 



and 



i,h 



if i / /i. 



Furthermore, we set C := Rp, where the vector p > has been defined 
in (22). Then Remark 5 after Theorem 2 together with the contraction prin- 
ciple [7] and simple algebraic calculations yields that, under the assumptions 
of Corollary 4, the sequence (dfcW(afc, s^^i, . . . , „))fcgN satisfies a sample 
path large deviation principle on Vq with normalizing sequence {bk)keN and 
good rate function 

inf /^™™'^(z). 

zeXlj^, w='I>(z-Cid) 

The contraction principle shows that this rate function is also the rate 
function in the large deviation principle for the sequence of stationary re- 
flected Brownian motions {^{'B/\/T)k — Cid))k^n with normalizing sequence 
{bk)keN- In particular, results concerning the rate function of reflected Brow- 
nian motion [2, 10, 15, 16, 18, 24] might be applied to obtain statements 
about the asymptotic behavior of tail probabilities of multiclass feedforward 
queueing networks in heavy traffic. 
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APPENDIX A: MODERATE DEVIATIONS FOR RENEWAL 

PROCESSES 

Here we consider for every G N a sequence (-'^fc,h)hez of independent and 
identically distributed strictly positive random variables with finite second 
moments. We derive moderate deviation principles for certain sequences of 
renewal processes constructed from these random variables in the topology 
II • II, from known moderate deviation principles for partial sums processes 
in the topology of uniform convergence on compacts ([33], Lemma 6.1). 

Lemma A.l. If there are values ^ > and S > such that 

lim E(Xfc,i)=e, 

fe— >oo 

lim Var(Xfc,i) = S^, 

fc— »oo ' 

and constants c > and 5 > such that, for every y £ [—6, 6] and k gN, 
Akiy) := \ogE{y{Xk,i - E{Xk,i))) < cy^ , 

then for every sequence {bk)keN of positive numbers satisfying (25) and (26), 
the centered and scaled sequence of linearly interpolated partial sums pro- 
cesses 

- xf-(fc-)-j^(Xfc,i)id(fc-) 

satisfies a sample path large deviation principle on Dq with good rate function 
jBrowny^2 normalizing sequence {bk)keN- 

Proof. We define for A; G N the piecewise constant partial sums (peps) 
process x^'^'''^ by setting X^'^''**(t) := xjj''**([tj) for t S M. Under the assump- 
tions of the lemma. Example 7.2 of [30] yields that the sequence of processes 

^X.l'^\k.)-EiXk,i)id{k-)- 



km 



Vbkk J keN 

satisfies a sample path large deviation principle on Dq with good rate func- 
tion /Brown ^'^2 ^^^^ normalizing sequence {bk)keN in the topology of uniform 
convergence on compacts. 

For every e > and r > 0, we have, for all A: G N with e^Jbkk > 

2supfcgN^(^fc,i)' 

/ X}^ fi 

P sup ' > e 

V-fcT-i<h</tr+i ybkk 

< {2kT + 2)P{Xk,i>eVhk) 

< {2kT + 2)P(exp((^(Xfc,i - EiXk,i)) - 6eVhk )) 

< (2fcT-F2)exp(c5^ -feVM?)- 
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Hence, for every e > 0, 

limsup — logP sup - — - — ^ — — ^ — ->e 

< lim sup — log P f sup ^ ^ 

/log(2fcr + 2) + c52 

<limsup "\ 1~ = — oo. 

fc^oo \ Ofc V J 

This implies by Theorem 4.2.13 in [7] that the sequence of processes 

'xf\k.)-E{Xk,i)id{k.y 



/bkk 



km 



also satisfies a sample path large deviation principle on Dq with good rate 
function /Brown^j-;2 ^^^^ normalizing sequence {bk)k€N in the topology of 
uniform convergence on compacts. 

Next, we strengthen this large deviation principle to the topology induced 
by the norm || • ||. According to Theorem 18 in [20], one can reach this goal 
by showing that, for every e,!^ > 0, there exists a r > such that 



lips^ 

>e \ <-v. 



1 / \-X.''^\kt) - E{Xk,i)kt\ 
lim sup — log F I sup — — 

fc^oo Ok \\t\>T Wbkk 

We choose e, > 0. A Chernoff type estimate gives, for every /i G N and 
7G[0,2c(5], 

p(j2(^k,m-E{Xk,l))>hj) 
\m=l / 

< Eexp(^^ (^k,m - EiXk,i)) - ^] 

\ m=l / 

= exp(/.A,(^)-^)<exp(-^). 
Similarly, one obtains, for every /i G N and 7 G [0, 2c5], 



logP (^k,ra - EiXk,l)) <-h^]< 



V m=l 



4c 
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Thus, we get, for every r G N and /c G N with e\Jbk/k < 2c6, 
\x'^'ikt)-E{Xk,i)kt\ 



sup- 

\t>T 



oo 

h=kr 
oo 

< 2exp 

h=kT 



t\/bkk 

h 



> e 



E{Xk,i)) 



m=l 



> he\ 



hbkS 



2exp(-T6fceV(4c)) 



( 4ck J 1 - expi^-bke"^ / {4ck)) ' 



This impHes that, for r G N with r > Acv /e^ 



hm sup 

fe— >oo ''fc 



1 / |X^'P^(fct)-^(Xfc,i)fct| 
— log P sup - — - — ^ — —^== — — - 

\t> 



< hm sup ( — 

fe— >oo V 



> e 



t^bkk 

log(l - expi-eHk/{4ck))) 



-V + hmsup 

fe— >oo 



4c bk 
log(4cA;/(e26fc)) 



bk 



-V. 



Here we have used the assumptions (25) and (26). A similar bound holds on 
the negative time interval which yields the envisaged strengthening through 
the principle of the largest term. □ 

We define the process by Xfc(t) := (x|jP'')-i(A:t) for t G R. Analogously 
to Lemma 2.4 in [32], we get the following: 

Corollary A. 2. Under the assumptions of Lemma A.l, the sequence 
of processes 

'^k-E{Xk,i)-^kid' 



^bkk 



feGN 



satisfies a sample path large deviation principle on Dq with good rate function 
^3jBrown^j-;2 normalizing sequence (bfc)^^^. 

For /c G N, we let be a random variable which is independent of the 
sequence (^fe,h)/iGZ and uniformly distributed on [0,1], and define the pro- 
cess x';j(t) := L(^fc''*^ ~ -^^fc,!)""*"] (kt). Clearly, we have |x'^ — x/d < 1, which 
yields by exponential equivalence [7] 

Corollary A. 3. Under the assumptions of Lemma A.l, the sequence 
of processes 



Xj_ 



E{X_ 



fc,l; 



^hk 



^kid\ 



km 
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satisfies a sample path large deviation principle on T>q with good rate function 
^3jBrowny'j-;2 normalizing sequence (6fc)fcGN- 

APPENDIX B: STATIONARITY OF NETWORK BEHAVIOR 

In this section we state independence and stationary increments condi- 
tions for random regular network primitives which imply the stationarity of 
the workload and queue length processes. 

We first investigate the behavior of a single queueing node under time 
shifts. For c G M and a function d G D, we let Ocd G P be the time shifted 
function defined by (Bcd)(i) := d{t + c) for every t G M. Furthermore, we set 
Scd := Gcd — d(c), which implies (Hcd)(0) = for every c G M and d G P. 
We use these shift mappings on product function spaces by applying them 
componentwise. For instance, if c G M", d G and i < n, we let the ith 
component of Hcd G be defined by (Hcd)j := Hc^dj. 

Lemma B.l. // (a, s) G x ij^^^ is a regular node primitives pair and 
t£R, 



Sa(t)(s"^) 


— ('='s-ioa(t)S) ) 


Hf(s~^oa) 


= (Ss-ioa(t)S)"^°(Sta), 


etV^(a,s) 


= V'^(Hja, Hs-ioa(t)S), 


GtW'^'^(a,s) 


= W^'^(Hta, H,-i„^(t)s), 


etQ«'^(a,s) 


= Q'^''^(Hja,Hs-ioa(t)s), 


et(Z^'^(a,s)-id) 


= Z'^''^(Hja, Hs-ioa(t)s) - 



and 

HtD^'^(a,s)=D^'^(Hta,H,-i„,(i)s) + Q^'^(Sta,H,-i„,(t)s)(0). 

Proof. Under the assumptions of the lemma, the first two lines in the 
display follow from the fact that, for every 6, c G M and d G linv, we have 
(Ocd — b)~^ = Gfe(d"-'-) — c. Next, one obtains the third and fourth line of 
the display by using the equation Oc(supd) — b = sup(0cd — 6), which is 
valid for every d G I'sup and c,b gM. For j GTC and c, t G M, we calculate 

-{ecQf^{a,s)){t) = nf^{a,s){t + c) - a,(t + c) 

= sup ^ji^) ~ + c) 

T<t+c, J2eeH ^7 ^°^<(^)^Efe-H s^^oa^ (t+c)-VW (a,s)(t+c) 

= sup (Hcaj)(T-c) 

^-^^*.E^gH("-(^r'°af))(r-c)<^^g^(H,(s7^oaf))(t)-e,VH(a,s)(t) 
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-(Hea,-)(t) 

'hi Q, 

= -Qj' (Sca,Hs-ioa(c)S)(t). 

This implies the fifth equation for components j GTC. The remaining equa- 
tions follow along similar transformations which are omitted. □ 

We use these shift properties to identify situations where the behavior of 
the feedforward network is stationary. 

Lemma B.2. We let (a, si, . . . , s„) el^ x Jj^^^^^ x • • • x X;^;^^" be a 
random regular network primitives tuple satisfying the following four condi- 
tions: 

1. The process a is independent of Si for every i <n. 

2. The processes Sij for i <n and j GTCiU Ci are independent. 

3. The process a has stationary increments, that is, the process Hja has the 
same distribution for every t G M. 

4. For every i<n and j G TliDCi, the process s~j has stationary increments 
for time shifts with values from the image of kj, which means that the 
distribution of the process Ea{s^j) is the same for every a E {aj(t) : t G 
R}CM. 

Then the processes Vi{a,si, . . . ,s„), . . . , V„(a, si, . . . , s„), Wi(a,si, . . . ,s„), 
W„(a,si,...,s„), Qi(a,si, . . . ,s„), Q„(a, si, . . . , s„), and Zi(a,si, 
. . . ,s„) - id, . . . , Z„(a, Si, . . . ,s„) — id are jointly stationary, that is, their 
common distribution is invariant under the shift Qt for every t £M. 

Proof. The conditions of the lemma together with the second equa- 
tion of Lemma B.l imply that the distribution of the pair (Hja, Hg-i^^j-j^si) 

does not depend on t G M. Hence, the joint stationarity of the processes 
Vi(a,si, . . . ,s„), Wi(a,si, . . . ,s„), Qi(a, si, . . . , s^) and Zi(a, si, . . . , s„) - 
id follows from the third to sixth equation of lemma B.l and the definitions 
of these processes. 

Setting di := Di(a, si, . . . , s„) = D^i'-^i (a, si), the last equation of Lemma B.l 
yields that the distribution of the triple (H^a, Hg-i^^^^^^si, Hjdi) does not de- 
pend on i G M. Clearly, the image of dij is a subset of the image of a^ for 
every j G A^. In view of the independence conditions 1 and 2, the process di 
is independent of the processes §2, . . . , s„. In particular, by Assumption 4, 
the distribution of the quadruple (Hia, Sg-i^^^^^^si, Htdi, H.-i^^^^^^S2) does 
not depend on t G M. Repeating the reasoning of the first paragraph simulta- 
neously for the first two queues, one therefore obtains the joint stationarity 
of the considered queueing processes of queues 1 and 2. 
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Iterating the scheme of the second paragraph through the remaining 
queues yields the statement of the lemma. □ 
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